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Radiative equilibrium surface temperatures and surface heating rates from a combined 
inviscid-boundary layer method are presented for the X-34 Reusable Launch Vehicle for 
several points along the hypersonic descent portion of its trajectory. Inviscid, perfect-gas 
solutions are generated with the Langley Aerothermodynamic Upwind Relaxation Algo- 
rithm (LAURA) and the Data-Parallel Lower-Upper Relaxation (DPLUR) code. Surface 
temperatures and heating rates are then computed using the Langley Approximate Three- 
Dimensional Convective Heating (LATCH) engineering code employing both laminar and 
turbulent flow models. The combined inviscid-boundary layer method provides accu- 
rate predictions of surface temperatures over most of the vehicle and requires much less 
computational effort than a Navier-Stokes code. This enables the generation of a more 
thorough aerothermal database which is necessary to design the thermal protection sys- 
tem and specify the vehicle’s flight limits. 


Introduction 

T HE X-34 1 6 is a reusable, sub-orbital test vehi- 
cle developed by Orbital Sciences Corporation 
(OSC) as part of NASA’s Reusable Launch Vehicle 
(RLV) Technology program. 7 Originally proposed as a 
partially reusable, two-stage vehicle designed to deliver 
1500 lb to low Earth orbit, the X-34’s current pur- 
pose is to provide a testbed for RLV technologies and 
to demonstrate RLV operations. These technologies 
include autonomous landing systems, low-cost avion- 
ics, advanced thermal protection systems (TPS), and 
composite airframe and propellant tanks. Relatively 
small in size, it is 58 ft long with a wing span of 28 
ft and a gross weight of approximately 45,000 lb. One 
of the program’s goals is to develop a vehicle capable 
of achieving Mach 8 flight and an altitude of 250,000 
ft. Part of NASA Langley Research Center’s role in 
the X-34 program is to assist OSC by performing both 
computational fluid dynamic (CFD) analysis on the 
vehicle as well as aerodynamic and aeroheating wind 
tunnel testing. One of Langley’s tasks is to provide 
OSC with CFD predictions of entry heating rates to 
be used for the TPS design. 
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The design of a vehicle’s TPS involves two areas. 
First, the maximum surface temperature along a tra- 
jectory defines which materials may be used for the 
TPS over different regions of the vehicle. Conversely, 
once the TPS materials have been chosen, flight limits 
are required to ensure the temperature limits of the 
materials are not violated for off-nominal trajectories. 
CFD assists in this area with solutions at or near the 
peak heating point (as defined by a stagnation heat- 
ing rate) of a trajectory. Second, the total heat load 
over the flight trajectory defines the thickness of the 
TPS materials. Full CFD is not appropriate. Long 
computer run times for individual Navier-Stokes (N-S) 
solutions prohibit running the many solutions required 
to define the heating along a trajectory. Instead, en- 
gineering codes such as MINIVER 8 are typically used 
to provide the complete time histories of surface tem- 
perature used to compute total heat load. 

Another approach to defining the surface heating 
along a trajectory is to use a combined inviscid- 
boundary layer method. Inviscid CFD solutions are 
less costly than N-S solutions; therefore, more points 
on a trajectory can be computed using the same com- 
puter resources. Also, engineering boundary layer 
methods such as the LATCH (Langley Approximate 
Three-dimensional Convective Heating) code 9 provide 
reasonably accurate heating rates over much of the 
vehicle (e.g., stagnation region, wind side, regions 
without flow separation) and run in minutes on desk- 
top workstations. However, care must be exercised at 
high-altitude, low Reynolds number conditions where 
the shock layer can not be divided into separate invis- 
cid and boundary layer regions due to viscous interac- 
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tions. At these conditions, an inviscid-boundary layer 
approach is inappropriate; and methods that treat the 
entire shock layer as viscous, such as viscous-shock- 
layer (VSL), 10 parabolized Navier-Stokes (PNS), 11 or 
N-S solvers are necessary. Inviscid-boundary layer 
methods are not meant to replace but complement 
benchmark CFD solutions in the area of TPS design. 
Although it may still be prohibitive to cover a tra- 
jectory in detail using this approach, heating rates 
computed at selected points on the trajectory can be 
used to calibrate the temperature time histories from 
an engineering method. 12 

This paper details the use of an inviscid-boundary 
layer method to compute the surface heating rates 
over the X-34 at several points along a representa- 
tive trajectory supplied by OSC. Inviscid, perfect-gas 
solutions are generated with the Langley Aerothermo- 
dynamic Upwind Relaxation Algorithm (LAL^RA) and 
the Data-Parallel Lower-Upper Relaxation (DPLUR) 
code. LATCH is used to compute the surface heating 
rates and radiative equilibrium temperatures. Com- 
parisons of the surface heating rates and temperatures 
are made with viscous, thin-layer N-S solutions from 
LAURA. 13 Maximum wind-side, lee-side, and wing 
leading edge temperatures are estimated as well. This 
w T ork is part of a collective effort at NASA Lang- 
ley to provide OSC w T ith the aerothermal information 
necessary to design the TPS for the X-34 vehicle. Ad- 
ditional data delivered to OSC by Langley include 
benchmark CFD solutions, 13 experimental aeroheat- 
ing, 14 and time histories of surface temperature. 12 

Geometry 

The full X-34 vehicle configuration (version 
X0001215) is shown in Fig. 1(a), and the geometry 
used for the inviscid solutions is shown in Fig. 1(b). 
Gaps in the elevons and between different TPS mate- 
rials are not modeled. Since LATCH requires a single 
block topology, the area aft of the wing trailing edge 
(including the vertical tail and body flap) is not in- 
cluded in the inviscid geometry. 

Trajectory 

The X 1004601 “no bounce” trajectory that is ana- 
lyzed is shown in Figs. 2-5. This “no bounce” trajec- 
tory is designed to eliminate the possibility of bouncing 
off of the atmosphere after reentry and is used as 
the reference trajectory for maximum wind-side and 
lee-side heating. The angle of attack varies from 25 
deg to 8 deg during the hypersonic descent portion of 
the trajectory. At the request of OSC, the following 
inviscid-boundary layer cases (2 ascent , 4 descent) are 
computed and are listed in Table 1. For this paper, 
the cases presented in detail are the M = 6.32, a = 
23 deg case (t = 330 sec) and the M = 6, a = 15.22 
deg case (t = 340 sec) because N-S solutions exist for 



a) Pull configuration used for viscous solutions. 



b) Partial configuration used for inviscid solutions. 
Fig. 1 X-34 configurations. 



these. Nonetheless, all inviscid-boundary layer solu- 
tions are used to compute the time histories of surface 
temperature used for the TPS design as detailed in 
Ref. 12. 

Computational Mesh 

The inviscid volume grid is obtained by truncating 
a viscous volume grid for the X-34 at the wing trail- 
ing edge. Because the viscous grid contains many grid 
points to resolve gradients, the inviscid grid is thinned 
and the points are redistributed to reduce unneces- 
sary clustering. The resulting inviscid volume grid is 
120 x 152 x 32 cells. The corresponding viscous vol- 
ume grid contains 64 cells between the body and grid 
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(deg) 


Table 1 Inviscid-boundary l ayer solutions. 

Thml Altitude Mach a S7i Inviscid N-S soln. Note 
(sec) (kft) No. (deg) (deg) code used available? 


145 183 6 9 0 LAURA 

152 196 6.83 11 0 LAURA 

* 330 118 6.32 23 0 DPLUR 

+10 LAURA 

* 340 112 6 15.22 0 DPLUR 

355 110 5.8 8 -10 LAURA 

578 86 3.6 6.46 -10 LAURA 



time ( sec ) 


Fig. 3 X1004601 trajectory (Mach number). 



time (sec) 


Fig. 4 X1004601 trajectory (angle of attack). 


No Mach 6 ascent 

No Max. heating on ascent 

Yes Max. heating 

Yes 

Yes Wind tunnel comparison 
No Min. a, max. heating 

No Reentry max, q, max, heating 



time (eec) 

Fig. 5 X1004601 trajectory (stagnation point 

heating). 

outer boundary. Grid resolution studies in Ref. 15 
indicate 32 cells is sufficient for inviscid calculations. 
Although DPLUR used this grid size for both of its 
solutions, LAURA’S multiblock capabilities allowed a 
coarser grid in the circumferential direction to be used 
in the nose region for its computations. This coarsen- 
ing of the grid speeds convergence of the solution in 
the stagnation region. Details of the grid generation 
process for the X-34 vehicle are given in Ref. 16. A 
sample, coarsened, flow-field grid is shown in Fig. 6. 

Flow-field Codes 

Inviscid solutions for the cases listed in Table 1 are 
generated with the CFD codes LAURA and DPLUR. 
To provide surface heating information to OSC in a 
timely fashion, two codes are used instead of one. 
Each code is tuned for a different computer architec- 
ture; LAURA for multitasking vector computers and 
DPLUR for massively parallel machines. Solutions 
for the database can be run concurrently on differ- 
ent systems which saves time. The inviscid flow fields 
serve as inputs to the LATCH engineering code which 
computes surface streamlines and both laminar and 
turbulent heating rates. Following are brief descrip- 
tions of the three methods. 
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Fig. 6 Inviscid flow-field grid (coarsened). 
LAURA 

LAURA (Langley Aerothermodynamic Upwind 
Relaxation Algorithm) is a finite-volume, shock- 
capturing algorithm for the steady-state solution of 
inviscid or viscous, hypersonic flows on rectangularly 
ordered, structured grids. LAURA has been used 
extensively to provide aerothermodynamic character- 
istics for a number of aerospace vehicles (e.g. AFE, 17 
HL-20, 18 Shuttle Orbiter, 19 Mars Pathfinder, 20 SSTO 
Access to Space 21 ) and is currently being used in 
the design and evaluation of the X-33 RLV. 22 The 
upwind-biased inviscid flux is constructed using Roe’s 
flux-difference-splitting 23 and Harten’s entropy fix 24 
with second-order corrections based on Yee’s sym- 
metric total-variation-diminishing scheme. 25 A point- 
implicit strategy is used which treats the variables at 
the cell center of interest implicitly at the advanced 
iteration level and uses the latest available data from 
neighboring cells. This results in an efficient, parallel 
implementation on multitasking vector computers. 26 
Gas chemistry options include perfect gas, equilibrium 
air, and air in chemical and thermal nonequilibrium. 
The algebraic turbulence models of Cebici-Smith 27 
and Baldwin-Lomax 28 are also available. More details 
of the algorithm can be found in Refs. 26, 29, and 30. 

DPLUR 

The DPLUR (Data-Parallel Lower-Upper Relax- 
ation) method 31,32 is based on the lower-upper sym- 
metric Gauss-Seidel method of Yoon and Jameson 33 
but has been modified for data-parallel computing. 
The Gauss-Seidel sweeps of the original method of 
Yoon and Jameson are replaced with a series of point 
Jacobi-like subiterations. This removes all data depen- 
dencies and yields a method that is almost perfectly 
parallel. Like LAURA, it is a finite- volume, shock- 
capturing algorithm for the steady-state solution of 


both inviscid and viscous flow fields on structured 
grids. Presently, there are options for perfect gas, equi- 
librium air, and 5-species nonequilibrium gas chem- 
istry. 

LATCH 

The engineering code LATCH (Langley Approxi- 
mate Three-Dimensional Convective Heating) 9 com- 
putes surface heating rates on three-dimensional (3-D) 
vehicles at angle of attack. The method is based 
on the axisymmetric analog for 3-D boundary layers 
and uses a generalized body fitted coordinate system. 
Boundary-layer edge conditions and the surface veloci- 
ties used to determine inviscid streamline direction are 
obtained from an inviscid flow-field solution. In this 
paper, inviscid solutions are supplied by both LAURA 
and DPLUR. Instead of solving the boundary-layer 
equations along streamlines, an approximate heating 
method developed by Zoby 34 that is valid for both 
laminar and turbulent heating is used. This method 
has been shown to produce accurate results for both 
wind tunnel and flight conditions 34-36 with only a frac- 
tion of the computational effort required by the full 
boundary-layer equations. 

Computational Resources 

The primary advantage to using an inviscid- 
boundary layer method over a N-S code is the reduc- 
tion in time needed to generate a solution. For the 
inviscid cases listed in Table 1, DPLUR requires 100 
node-hours per solution on an IBM SP-2 and LAURA 
requires about 25 hours per solution on a CRAY YMP. 
LATCH boundary layer solutions containing surface 
temperatures and heating rates are then obtained in 
about 5 minutes each on an SGI R10000 worksta- 
tion. Conversely, the viscous LAURA solutions require 
approximately 300 hours each on a CRAY C-90 to 
reach convergence. Although the inviscid-boundary 
layer approach still uses a considerable amount of com- 
puter time and resources compared to pure engineering 
methods, it offers a significant savings over viscous N-S 
CFD codes. 


Results 

Surface temperature contours and heating rates are 
examined for the X-34 at M = 6.32, a =23 deg 
(t = 330 sec) and M = 6, a =15.22 deg (t = 
340 sec). Results from a combined inviscid-boundary 
layer method (LAUR A-LATCH and DPLUR-LATCH) 
are compared with viscous solutions from LAURA to 
assess the accuracy of the inviscid-boundary layer ap- 
proach. Both laminar and turbulent solutions are com- 
puted although OSC only requested turbulent heating 
data. All solutions assume a perfect gas and com- 
pute radiative equilibrium wall temperatures based on 
an emissivity of 0.8. The turbulent, viscous solutions 
from LAURA are computed using the Baldwin-Lomax 
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Fig. 7 Laminar wind-side temperatures at t = 
330 sec. 



Fig. 8 Turbulent wind-side temperatures at t = 
330 sec. 


Turbulent WindskJe 



Fig. 9 Turbulent wind-side temperatures at t = 
330 sec (with elevons deflected). 

algebraic turbulence model. Details of the LAURA 
viscous solutions are found in Ref. 13. 

Wind-side temperature contours are shown in 
Figs. 7-10 for the two flight conditions. Contour lev- 
els are plotted in 100 deg F increments over a range 
of temperatures from 300 2000 deg F. The upper half 
of each figure depicts the inviscid-boundary layer re- 
sults (either LAURA-LATCH or DPLUR-LATCH), 
and the lower half shows the surface temperatures 
from a LAURA viscous solution. To help correlate 



Fig. 10 Turbulent wind- side temperatures at t = 
340 sec. 


Table 2 Multi-use temperature limits of TPS blan- 
kets. 


Material 

Max. Temperature (deg F) 

High Temp. AFRSI 

2000 

AFRSI 

1500 

FRSI 

700 


the predicted surface temperatures with TPS materi- 
als, Table 2 lists the multi-use capability of the TPS 
blankets used over much of the vehicle. In Fig. 7, 
the laminar temperatures from DPLUR-LATCH agree 
quite well (i.e. within 100 deg F) with the tempera- 
tures from LAURA over much of the lower surface. 
Both solutions predict temperatures of 800 deg F near 
the centerline, 1100 deg F near the outer edge of the 
strake, and 900-1000 deg F near the middle of the 
wing. This region of higher temperatures extending 
across the wing from the leading to the trailing edge 
is caused by the wing-bow shock interaction and is 
predicted by both methods. However, the magnitude 
of the temperatures appears to be slightly lower for 
DPLUR-LATCH. Similar comparisons are seen for the 
turbulent results shown in Figs. 8-10, albeit the overall 
temperature levels are 300-600 deg F higher than the 
laminar temperatures. The higher temperatures force 
the use of the High-Temperature AFRSI blankets over 
much of the windward surface. Figure 9 shows the 
temperature contours for the X-34 with +10 deg de- 
flected elevons. Both LAURA-LATCH and LAURA 
show a pocket of higher temperatures 1700 deg F) 
on the elevon surface. In addition, Figs. 8 and 9 offer 
a comparison between DPLUR-LATCH and LAURA- 
LATCH solutions at the same conditions (except for 
the elevon deflection). The temperature contours for 
the two solutions are similar over much of the lower 
surface except near the forward portion of the vehicle 
where DPLUR-LATCH predicts a lower temperature 
awny from the centerline. As previously stated, a 
coarser grid was used for the inviscid LAURA solution 
in the circumferential direction in this region which 
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Laminar Laaside 



Fig. 11 Laminar lee-side temperatures at t — 
330 sec. 


Turbulent Laaatde 



Fig. 12 Turbulent lee-side temperatures at t = 
330 sec. 

Turbulent Laaside 



Fig. 13 Turbulent lee-side temperatures at f = 
340 sec. 

may account for some of the differences. Figure 10 
shows the wind-side temperature contours for t = 340 
sec. As expected, the lower angle of attack (a =15,22 
deg) results in generally lower temperatures on the 
wind side of the vehicle. For example, the radiative 
equilibrium temperatures over most of the wing are 
around 200 deg F lower than for the t = 330 sec case. 

Leeside temperature contours are shown in Figs. 11- 
13. Temperature contours for the turbulent, deflected 
elevon case at t = 330 sec are not presented since the 


Laminar Side 
T*, deg F 



LAURA (Viscous) 


a -23 deg 
S„ = 0 deg 

Fig. 14 Laminar side temperatures at t — 330 sec. 


Turbulent Side 
T* deg F 



LAURA (Viscous) 


a = 23 deg 
5* = 0 deg 

Fig. 15 Turbulent side temperatures at f = 
330 sec. 

lee-side temperatures on the elevon are very low (< 
300 deg F). Somewhat unexpectedly for an inviscid- 
boundary layer method, the lee-side temperature con- 
tours from DPLUR-LATCH agree quite well with the 
viscous LAURA solution at these conditions. The 
same general patterns in temperature are seen near 
the forward portion of the vehicle as well as on the 
wing. These lee-side predictions impact the TPS de- 
sign because the temperature levels vary around the 
limit of the FRSI blankets (700 deg F). 

Temperature contours on the side of the vehicle are 
shown in Figs. 14-16. Again the contour patterns 
from DPLUR-LATCH and the viscous LAURA solu- 
tions are very similar. Therefore, as seen in Figs. 7- 16, 
the inviscid-boundary layer technique predicts radia- 
tive equilibrium wall temperatures that compare fa- 
vorably with temperatures from a N-S solver. 

To illustrate the differences between the inviscid- 
boundary layer and N-S solutions more clearly, surface 
heating rates are examined along several cut planes 
including centerline, wing leading edge, and cross- 
sectional cuts. Surface heating rates are more sen- 
sitive than radiative equilibrium wall temperatures 


r 


* 
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Turbulent Side 
T* degF 



M _ 6 LAURA (Viscous) 

as 15.22 deg 
5- = 0deg 


Fig. 16 Turbulent side temperatures at < = 
340 sec. 



Fig. 17 Cut plane locations. 

(q w ~ T£) and should provide more insight into the 
comparison between the methods. The locations of 
cut planes are shown in Fig. 17. 

Lateral surface heating distributions are given in 
Figs. 18 20 at four axial stations. The laminar heating 
rates in Fig. 18 are examined first, because there is no 
influence of different turbulence models. The laminar 
cuts in Fig. 18 show, in general, an overall good agree- 
ment between the DPLUR-LATCH, LAURA-LATCH, 
and the viscous LAURA solutions. However, the 
LATCH results tend to underestimate the maximum 
heating rates from LAURA by approximately 15-20 
percent at each axial station. Because LATCH com- 
putes heating rates along inviscid surface streamlines, 
it has difficulties in regions of high curvature such 
as near the wing leading edge. The surface stream- 
lines cannot account for the large three-dimensional 
effects that are present. In particular, laminar heat- 
ing rates from LATCH in high-curvature regions tend 
to be lower than those predicted by N-S solvers. 9 
This may help explain the heating rate comparisons 
in Figs. 18(c) and 18(d) and to a lesser degree in 
Figs. 18(a) and 18(b) where the streamlines wrap from 
the lower surface around to the side of the vehicle. 

Differences in the heating pattern near the wind- 


ward centerline are noted in Figs. 18(a) and 18(b) 
between LATCH and LAURA. The heating decreases 
- as the centerline is approached. A possible explanation 
is suggested by the DPLUR-LATCH and LAURA- 
LATCH results in Fig. 18(b). LAURA-LATCH, which 
uses a coarser circumferential grid in this region, does 
not predict the dip in heating as well as DPLUR- 
LATCH, which uses the full grid. The streamline 
directions on the flat lower surface of the vehicle 
are sensitive to the grid resolution. Also seen in 
Figs. 18(b) and 18(c) is the increased heating predicted 
by LAURA near the leeward centerline. Crossflow sep- 
aration is present around x = 150 in on the lee side of 
the vehicle, 13 and this is reflected in the higher heat- 
ing rates. Being an engineering code, LATCH cannot 
account for this. 

The turbulent heating rates in Figs. 19 and 20 are 
examined next. In general, good agreement between 
LATCH and LAURA is seen in the cross-sectional 
cuts, especially on the wing at x = 300 in and x 
— 450 in. However, unlike the laminar results, the 
turbulent peak heating rates from LATCH tend to 
be higher than those from LAURA. This fact can be 
attributed to the inherent differences between the engi- 
neering turbulent boundary-layer equations in LATCH 
and the algebraic model employed by LAURA. For the 
turbulent cases, the heating patterns near the wind- 
ward centerline in Figs. 19(b) and 20(b) differ between 
LATCH and LAURA. Vortices on the wind side of the 
vehicle in the boundary layer have been observed in 
the flow fields predicted by LAURA that might explain 
these differences. 13 However, to keep the comparisons 
in perspective, differences of 15-20 percent in heating 
rates with corresponding differences of 4 5 percent in 
temperature are adequate for design work. 

Windward centerline distributions for the two cases 
(t = 330 sec and t = 340 sec) are presented in Fig. 21. 
The turbulent heating rates from LATCH (with both 
LAURA and DPLUR) are approximately 15 percent 
lower than the heating rates from LALTRA for most of 
the vehicle (x > 100 in) for both cases. The laminar 
heating rates in Fig. 21(a) from LATCH are 25 percent 
higher than the LAURA results for the forward half 
of the vehicle (x < 300 in). The agreement is much 
better downstream. The heating rates from DPLUR- 
LATCH are closer to the LAURA heating rates due 
to the denser circumferential grid resolution used by 
DPLUR. 

Wing leading edge heating distributions are pre- 
sented in Figs. 22-24. The leading edge is defined as 
the outermost point of the wing (z = z maXl see Fig. 17) 
and does not necessarily represent the highest heat- 
ing rates or temperatures on the wing. In Figure 22, 
the laminar heating rates from LAURA are approxi- 
mately 15 percent higher than the invisc.id-boundary 
layer '(DPLUR-LATCH and LAURA-LATCH) solu- 
tions along the wing leading edge. This corresponds to 
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Fig. 18 Lateral laminar heating distributions at t = 330 sec. 


peak temperatures (not shown) along the leading edge 
of 1800 deg F for LAURA and 1725 deg F for LATCH. 
As shown previously, it is the inviscid-boundary layer 
solutions that predict the higher turbulent heating 
rates as shown in Figs. 23 and 24. Turbulent peak 
temperatures are 2035 deg F (t = 330 sec) and 19G5 
deg F (t = 340 sec) from LAURA and 2110 deg F 
(t = 330 sec) and 2020 deg F (t = 340 sec) from 
LATCH. SIRCA tiles are used for the wing leading 
edge, because these temperatures exceed the High- 
Temperature AFRSI limit of 2000 deg F. Overall, the 
inviscid-boundary layer method predicts reasonably 
good (within 15 percent) surface heating rates and 
radiative equilibrium surface temperatures along the 
wing leading edge, especially considering the strong 
bow'-wing shock interaction. 


Wing heating distributions at z = 100 in. are pre- 
sented in Figs. 25-27. Figures 25 and 26 show 
the effect of elevon deflection on the heating rates. 
The inviscid-boundary layer results agree well with the 
viscous LAURA solutions over much of the wing in- 
cluding the elevon. Heating rates from LATCH are 
approximately 20 percent lower than LAURA for the 
laminar case presented in Fig. 25 but are within 10 
percent for the turbulent cases. Although not shown, 
the deflected elevon results in increases of about 190 
deg F for laminar flow and 240 deg F for turbulent 
flow over the surface temperatures of the undeflected 
elevon. Both LATCH and LAURA predict similar 
jumps in temperature. 


8 OF 12 


American Institute of Aeronautics and Astronautics Paper 98 0880 




a) x = 50 in. 


b) x = 150 in. 



z (in) 


c) x = 300 in. 



Fig. 19 Lateral turbulent heating distributions at t = 330 sec. 


Concluding Remarks 

A combined inviscid-boundary layer method 
(LAURA-LATCH, DPLUR-LATCH) has been used 
to predict the surface heating rates and radiative 
equilibrium wall temperatures for the X-34 vehicle 
along a reference trajectory. This information has 
been delivered to OSC as part of a collective effort 
by NASA Langley to aid the TPS design. Wall 
temperature patterns from the engineering boundary 
layer code LATCH are similar to the wall tem- 
peratures from the N-S solver LAURA over much 
of the vehicle at two flight conditions. Increased 
temperatures along the wing due to the wing-bow 
shock interaction and on the deflected elevon are 
correctly predicted by the inviscid-boundary layer 
technique. LATCH predicts surface heating rates 


that are generally within 20 percent of values from 
a viscous LAURA solution. The observed agreement 
between LATCH and LAURA is somewhat better for 
turbulent flows. The turbulent radiative equilibrium 
surface temperatures are 300 600 deg F higher than 
the corresponding laminar temperatures at the same 
conditions. The inviscid-boundary layer method 
(DPLUR-LATCH and LAURA-LATCH) also uses 
much less computer time than the N-S solver LAURA. 
Consequently, many more solutions can be computed 
along a vehicle’s trajectory with the same computa- 
tional resources. This ability to generate a reasonably 
accurate aerothermal database for a vehicle makes 
inviscid-boundary layer methods excellent design 
tools. 
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c) X = 300 in. d) x = 450 in. 

Fig. 20 Lateral turbulent heating distributions at t = 340 sec. 
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